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The detection of a stochastic background of gravitational waves could significantly impact our un- 
derstanding of the physical processes that shaped the early Universe. The challenge lies in separating 
the cosmological signal from other stochastic processes such as instrument noise and astrophysical 
foregrounds. One approach is to build two or more detectors and cross correlate their output, 
thereby enhancing the common gravitational wave signal relative to the uncorrelated instrument 
noise. When only one detector is available, as will likely be the case with the Laser Interferome- 
ter Space Antenna (LISA), alternative analysis techniques must be developed. Here we show that 
models of the noise and signal transfer functions can be used to tease apart the gravitational and 
instrument noise contributions. We discuss the role of gravitational wave insensitive "null channels" 
formed from particular combinations of the time delay interferometry, and derive a new combination 
that maintains this insensitivity for unequal arm length detectors. We show that, in the absence of 
astrophysical foregrounds, LISA could detect signals with energy densities as low as fi gw = 6 x 10 -13 
with just one month of data. We describe an end-to-end Bayesian analysis pipeline that is able to 
search for, characterize and assign confidence levels for the detection of a stochastic gravitational 
wave background, and demonstrate the effectiveness of this approach using simulated data from the 
third round of Mock LISA Data Challenges. 
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I. INTRODUCTION 

One of the most exciting prospects for the Laser Inter- 
ferometer Space Antenna (LISA) observatory is the pos- 
sibility of detecting a stochastic gravitational wave back- 
ground from the early Universe. Just as studies of the 
cosmic microwave background have revolutionized our 
understanding of cosmology [l[ , the detection of a grav- 
itational wave background would provide unique insight 
into the processes that shaped the early universe 043 • 

The challenge lies in separating a stochastic back- 
ground from instrument noise and foreground sources. 
When multiple independent detectors are available, as 
is the case with the ground based interferometers, the 
signal from one interferometer can be used as a (noise 
corrupted) template for a second interferometer. The 
common gravitational wave signal will combine coher- 
ently, while the contributions from instrument noise will 
average to zero 

With a single detector, such as LISA, the cross cor- 
relation technique will not be available. The key to de- 
tecting a stochastic background with LISA is the devel- 
opment of a detailed understanding of how instrument 
noise and gravitational wave signals manifest in the vari- 
ous interferometry channels. Here we show how a model 
for the six cross-spectra of the Time Delay Interferom- 
etry (TDI) [O, [H| channels can be used to tease apart 
a stochastic background from instrument noise. Key to 
this approach are the very different geometrical transfer 
functions that affect the instrument noise and a stochas- 
tic gravitational wave background. While the instrument 
noise spectra could in principle conspire to mimic the ge- 
ometrical transfer functions of the signal, or vice-versa, 
this is unlikely to occur in practice as the signal and noise 
models have very distinct and complicated transfer func- 



tions. The robustness of our approach can be enhanced 
by developing informative priors for the noise and sig- 
nal spectra. This can be done through a combination 
of preflight testing, on-orbit commissioning studies, and 
theoretical modeling. We use the TDI cross-spectra to 
compute a likelihood function, which is multiplied by the 
priors to yield the posterior distributions for the param- 
eters describing the instrument noise and cosmological 
background. 

Models for the TDI cross-spectra have previously been 
considered in the context of LISA instrument noise de- 
termination flil ]. and the insensitivity of certain TDI 
variables to gravitational wave signals have been put 
forward as a technique for discriminating between a 
stochastic background and instrument noise [ill, EH EH ■ 
Our approach extends the noise characterization study 
of Sylvestre and Tinto [14| to include the signal cross- 



spectra, and improves upon the simple estimator used by 
Hogan and Bender [l5| by using an optimal combination 
of all six cross-spectra. Our approach is able to detect 
stochastic signals buried well below the instrument noise. 

In the remainder of the paper we describe our noise and 
signal models for both equal arm and unequal arm LISA 
interferometers, and develop an end-to-end Bayesian 
analysis pipeline that is able to search for, character- 
ize and assign confidence levels for the detection of a 
stochastic gravitational wave background. We apply our 
approach to simulated data from the third round of Mock 
LISA Data Challenge, and show that we are able to ac- 
curately recover the stochastic signal and independently 
measure the position and acceleration noise levels in each 
arm of the interferometer. We carry out our analysis us- 
ing a noise orthogonal A, E, T set of TDI channels. In 
this basis, and for equal arm-lengths, the T channel is in- 
sensitive to gravitational wave signals at low frequencies. 
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For the first time, we show that a modified signal insen- 
sitive T channel can be found for unequal arm-lengths. 

We close with a discussion of how our approach can be 
extended to account for the galactic foreground. While 
several tens of thousands of galactic sources can be in- 
dividually resolved and regressed from the data, the re- 
mainder will form a confusion noise for LISA [l6[ . We dis- 
cuss how the anisotropy of this foreground signal, along 
with prior information about the spectral amplitude pro- 
vided by the resolved systems, can be used to separate the 
galactic foreground from an isotropic background, but 
defer a detailed analysis to a future publication. 



II. EQUAL ARM LISA 
A. Noise Model 

The LISA constellation will be composed of three satel- 
lites in an approximately equilateral triangle configura- 
tion. Laser beams will be transmitted between each pair 
of satellites, and interferometry signals will be formed 
using these beams. Each of the three satellites will have 
two proof masses, one for each of the two incoming laser 
beams. We are looking to detect variations in the light 
travel time, or equivalently, the distance, between the 
proof masses along each arm of the interferometer. Noise 
enters the measurement when the proof masses move in 
response to local disturbances, and in the process of mea- 
suring the phase of the laser light. The various LISA 
noise sources are discussed in several references 

As is commonly done, we group all the noise sources 
into two categories, position and acceleration. Each proof 
mass will have a position and an acceleration noise asso- 
ciated with it, making a total of six position and six accel- 
eration noise levels. These twelve noise levels will be the 
parameters in our noise model. We assume that we un- 
derstand the noise spectra of each of the instrument com- 
ponents, through some combination of pre-launch test- 
ing, on-board instrument characterization and theoreti- 
cal modeling. In reality it may be necessary to include 
additional parameters in the noise model to account for 
uncertainties in the spectral shape of the individual noise 
spectra. However, if the priors on these additional pa- 
rameters are narrowly peaked, there will be little impact 
on our ability to detect a stochastic background signal. 

For this study we adopted the model used for the Mock 
LISA Data Challenges (MLDC) [20j]. The position noise 
affecting each proof mass is assumed to be white, with a 
nominal spectral density of 



of 



S a (f) = 9 x 10" 



(2) 

The precise level of each contribution is to be determined 
from the data. A more realistic model for the noise con- 
tributions would include a parameterized model for the 
frequency dependence, with the model parameters to be 
inferred from the data. Allowing for this additional free- 
dom would weaken the bounds that can be placed on the 
contribution from the stochastic background. 

We derive here the transfer functions describing how 
these twelve noise levels enter the data stream. We start 
by writing down the phase output $y(£), for the link 
connecting spacecraft i and j: 



*y(t) =Ci(t- Li-) - Cj{t) + iM*) 



Lij)) 



(3) 
(4) 



Here the Ci are the laser phase noises, ipy is the gravi- 
tational wave strain, and nF- and nfj denote the position 
and acceleration noise. The laser phase noise is canceled 
by using the Time Delay Interferometry(TDI) variables 
developed by Armstrong, Estabrook, and Tintoflll [l^|. 
The gravitational wave strain is assumed to be uncorre- 
lated with the noise and will have an expectation value 
of zero when multiplied by anything other than itself. 

A Michelson signal can be formed at any of the three 
vertices by combining the phase at that detector with the 
time delayed signal from the two detectors at the ends of 
the two adjacent arms: 

Afi(t) = $ia(t-Iia)+*ai(*)-*i2(*-ii3)-*3i(t) (5) 

The TDI channels which cancel laser phase noise for an 
equal arm LISA (Lij = L = 5 • 10 6 km) arc formed by 
subtracting a time delayed Michelson signal as follows: 



X(t) = Mi(t)-Mi{t-2L) 



(6) 



Moving to the frequency domain, the signal at vertex 1 
can be written as 
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(7) 



S p (f) = 4 x 10" 42 Hz" 



(1) 



The acceleration noise is taken as white above 0.1 mHz, 
with a red component below this frequency. Integrated 
to give an effective position noise, the proof mass distur- 
bances on each test mass have a nominal spectral density 



Here ///, = c/{2nL). The other TDI channels, Y and Z, 
are given by permuting indices in the expression for the 
X channel. Alternative combinations of the phase mea- 
surements can be used to derive noise orthogonal TDI 
variables [l3| . Here we use a different set of noise orthog- 
onal channels formed from combinations of the X.Y, and 
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Z channels: 



Equal-Arm LISA Data vs. Model 



A = -(2X -Y-Z) 
T = \(X + Y + Z) 



(8) 



We calculate the six cross-spectral densities for these 
channels in the appendix. As an example, we quote 
here the position and acceleration noise contributions to 

(AA*y. 



sm 2 f L 
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+5(S 21 + Sf 2 + 5f 3 + S 31 ) 

+2(5f 3 + 5f 2 )l (9) 
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where £y(/) 



j(/)n* (/)). Our model for the noise 



depends on the twelve noise spectral density levels, which 
we assume fully describe the LISA instrument noise. 
Note that the sums of the noise contributions in each 
arm of the interferometer, such as (S 23 + S 32 ), that are 
strongly constrained. The differences arc very weakly 
constrained by the model. Figure 1 compares our model 
spectra to simulated LISA noise spectra from the MLDC 
training data. 



B. Stochastic Background Model 

In this section we re-derive the LISA response to a 
stochastic gravitational wave background. We later ex- 
tend this calculation to take into account the effect of 




FIG. 1: Our model for the noise in the A and T channels com- 
pared to smoothed spectra formed from the MLDC training 
data. 



unequal arm-lengths. We begin by expanding the gravi- 
tational wave background in plane waves: 



(11) 

Here the e^- are the components of the polarization tensor 
and P sums over the two polarizations. The polarization 
tensors are formed by using the basis vectors u and v and 
the sky location vector Cl, 

u = cos 6 cos 4>x + cos 9 sin <jyy — sin 6z 
v = sin <j)x — cos <py 

ft = sin 6 cos <j>x + sin 9 sin 4>y + cos 62 (12) 
The polarization tensors are formed as 

e + (A,V) = e + (f7)cos(2i/>) - e x (£l) sin(2^) (13) 

and 

e x (Cl, V>) = e + (0) sin(2?/>) + e x (A) cos(2V>) (14) 

where tp is the polarization angle and e + and e x are given 
by: 



and 



u <X> u — v (K> v 



u <x> v + v <x> u . 



(15) 



(16) 



The response registered in an interfcromctry channel 
can be written as 



S i (t)=T> i (n,f):h(f,x) 



(17) 



where 



(hi ®hi)T(ra ■ k,f) 
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and 

r = 



Equal-Arm LISA Sensitivity Curv 
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The signal cross spectra are given by 

(S i (f),S j (f)) = S h (f)R ij (f) 



where 



and 
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(18) 



(19) 



(20) 



Here £l gw (f) is the energy density in gravitational waves 
per logarithmic frequency interval, scaled by the clo- 
sure density. We assumed a simple power law behavior: 
^gw(f) = (//lmHz) n fl 9W , where £l gw denotes the level 
at lmHz and n is the spectral index. 

The beam pattern functions are given by 

Ff(Cl,f)=D i (QJ):e p (Cl). (21) 

In general, the integral in (|19[) must be performed nu- 
merically, though we can develop analytic expressions in 
the low frequency limit: 



Raa = Ree = 4 sin 2 4" 
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The A,E, and T channels were created to be noise or- 
thogonal, but they also happen to be signal orthogonal 
in the equal arm case. All of the cross terms in the re- 
sponse function Rae, Rat etc. are zero in the equal 
arm-length limit. Sensitivity curves for the various chan- 
nels are generated by plotting 



h K = 



Skk (/) 
RkkU) 



(24) 



where Skk and Rkk are the noise and signal spectral 
densities in the K channel. Figure 2 shows the sensi- 
tivity curves for the A, E and T channels along with 
a scale invariant gravitational wave background with 



a, 
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We see that T channel is insensitive to 



the gravitational wave background for / < /* 



A,E Channel 
T Channel 
!X,„ = 1.0'10/ ra 



FIG. 2: The sensitivity curve for the A, E and T channels, 
showing the insensitivity of the T channel to a gravitational 
wave signal. 



III. SEARCH AND CHARACTERIZATION 

For Gaussian signals and noise, the likelihood of mea- 
suring cross-spectra < XiXj > is given by: 



P (x\x) = n 



i 



(2tt) n / 2 \C\ 



exp (XiC^Xj) 



(25) 



where C is the noise correlation matrix, Xi — {A, E, T}, 
N is the number of samples in each channel and x — > 
(S 1 ", S?, Qg W , n) denotes the parameters in our model. 
The noise correlation matrix is given by: 



f(AA) (AE) (AT)\ 
Ca = (EA) (EE) (ET) 
\ (TA) (TE) (TT) ) 



(26) 



Combined with a prior for the model parameters, p(x), 
we are able to generate samples from the posterior dis- 
tributions function p(x\X) = p(X.\x)p(x)/p(K) using 
a Parallel Tempered Markov Chain Monte Carlo algo- 
rithm (2~j| . Lacking a detailed instrument model or rele- 
vant experimental data, we choose to use uniform priors 
for the 12 instrument noise levels, allowing a factor of ten 
variation above and below the nominal levels. The en- 
ergy density was taken to be uniform in ]n(Cl gw ) across 
the range [—30,-24.5]. We considered two models for 
the spectral slope, cither assuming a scale invariant back- 
ground and fixing n = 0, or allowing the spectral index 
n to be uniform in the range [—1,1]. 

For the proposal distribution we used a mixture of uni- 
form draws from the full prior range, and draws from 
a multivariate Gaussian distribution computed from the 
Fisher Information Matrix 211. Correlations between the 



parameters, and the frequency dependence of the spectra 
complicate the computation of the Fisher Matrix, but the 
basic idea can be understood by considering zero mean 
white noise with variance a. The relevant question is, 
how well can the noise level a 2 be determined from N 
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noise samples? The likelihood of observing the data {x{\ 
is 



C 



1 



{2n) N / 2 a 



N 



exp 



2a 2 



(27) 



which yields a maximum likelihood estimate for the noise 
level of ct 2 = J2i x i l N - Tlie Fisher Matrix has a single 
element: 



d 2 ln£ , N 

J- era o o ML ^ ^ n 



(28) 



Thus, the error in the estimated noise level is Act 2 = 
2a Act = ct 2 / ' y/N/2. We see that the fractional er- 
ror in the noise level estimate scales as the square root 
of the number of data points. The same is true for 
the more complicated colored spectra in our LISA noise 
model. Noting that the acceleration noise dominates be- 
low ~ lmHz, while the position noise dominates above 
~ 3mHz, it follows that the effective number of samples 
available to constrain the acceleration noise is of order 
N a r~j lmHz x T bs- The position noise is far better con- 
strained, with order N p ~ /jv x T b s samples, where /at 
is the Nyquist frequency of the data. As mentioned ear- 
lier, only the sum of the instrument noise levels in each 
arm are strongly constrained, so the Fisher matrix ap- 
proach leads to very large jumps being proposed in the 
noise level differences in each arm. To maintain a good 
acceptance rate, we capped the variance in the weakly 
constrained directions to be ten times the variance in the 
well constrained directions. 



IV. MOCK LISA DATA CHALLENGE 

We are able to test our analysis technique on simu- 
lated data from the third round of the Mock LISA Data 
Challenge. Specifically, Challenge 3.5 provides month 
long training and blind data sets of 2 20 samples with 2 
second sampling. An isotropic gravitational wave back- 
ground was injected with a level that is approximately 10 
times the nominal noise levels at lmHz, corresponding to 
n GW = 8.95 x 10" 12 - 1.66 x 10" 11 for a Hubble con- 
stant of Hq = 70km/s/Mpc. The noise levels were drawn 
from a range within ±20% of the nominal values. Note 
that we used prior ranges far wider range than this as we 
wanted to test our approach in a more realistic setting. 



A. Training Data Results 

We found that the signal transfer functions in the 
training data did not match our analytic model. We 
traced the problem to the time domain filters that were 
used to generate the data sets, which introduced addi- 
tional transfer functions in the frequency domain. Since 
the analytic form for these transfer functions have not 
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FIG. 3: Histograms showing the posterior distribution func- 
tions for the position noise levels, scaled by the nominal level. 
One the left are the sums along each arm, and on the right are 
the differences. The vertical lines denote the injected values. 
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FIG. 4: Histograms showing the posterior distribution func- 
tions for the position noise levels, scaled by the nominal level. 
One the left are the sums along each arm, and on the right are 
the differences. The vertical lines denote the injected values. 
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FIG. 5: The posterior distribution for the gravitational wave 
background level (scaled up by 10 13 ). The vertical line de- 
notes the injected values. 
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FIG. 6: The posterior distribution for the gravitational wave 
background slope. 



been published, we used the training data to estimate 
the transfer functions and update our signal model. 

We used the Bayesian analysis algorithm described 
above and recovered the PDF for the noise parameters 
and stochastic background energy density show in Fig- 
ures 3-6. As mentioned earlier, only the sums of the noise 
contributions in each arm are constrained (the accelera- 
tion and position noise contributions can be separated 
though as they have very different transfer functions). 
Figure 3 provides an example of this by showing that 
the sums 5^ + Sfj are well constrained, while the differ- 
ences — Sfj are poorly constrained. Also note that 
the position noise levels are far better determined than 
the acceleration noise levels, as expected from our Fisher 
Matrix analysis. 
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FIG. 7: Bayes factor showing the detectability versus back- 
ground level. The line shows the average values of the scat- 
tered points. 



B. Detection Limits 

Having established that our algorithm can faithfully re- 
cover a stochastic background level of il gw ~ If) -11 with 
one month of data, we next turned to the problem of de- 
termining the LISA detection limit for an approximately 
scale invariant stochastic gravitational wave background. 
To do this we generated a new set of simulated data sets 
by re-scaling the gravitational wave contribution to the 
MLDC training data, and used Bayesian model selection 
to compared the evidence for two models, Mq - the data 
is described by instrument noise alone, and Mi - the data 
is described by instrument noise and a stochastic grav- 
itational wave background. Following the approach de- 
scribed in Ref. [2l|, we computed the evidence for each 
model using thermodynamic integration (22j . Plots of 
the Bayes factor, or evidence ratio, as a function of Q gw 
are shown in Figure 7. We performed multiple runs with 
different random number seeds as a way to estimate the 
numerical error in our Bayes factors. Our detection con- 
fidence becomes very strong (a Bayes factor of 30) for 
a background level of f2 gw = 6 x 10~ 13 with one month 
of data. With one year of data the limit improves to 
figw = 1.7 x 10~ 13 . 

It is interesting to note that the detection limit does 
not change if we include the spectral slope of the back- 
ground as a model parameter. At first this seems a little 
surprising, as the simulated data has a spectral slope of 
n = 0, so we would expect the simpler model with n = 
to be favored over the more complicated model with n as 
a free parameter. Further investigation revealed that the 
more complicated model was able to provide a slightly 
better fit to the data, and that this was enough to com- 
pensate for the additional complexity of the model in the 
calculation of the model evidence. Our hypothesis is that 
our signal model is imperfect because of the need to use 
numerical fits to the transfer functions introduce by the 
simulation software, and that the freedom to adjust the 
spectral slope of the background is able to compensate 
for this imperfection. 



C. The role of null channels and 4-link operation 

Previous analyses of stochastic background detection 
with LISA have emphasized the importance of the gravi- 
tational wave insensitive null channel that can be formed 
from the Sagnac or Michelson intcrfcromctry channels 
(to be precise these channels arc only null in the zero fre- 
quency limit, at non-zero frequencies they do respond 
weakly to gravitational signals because of finite arm- 
length effects). Hogan and Bender showed how this null 
channel could be used to construct a statistic that mea- 
sures the amplitude of the stochastic background. The 
importance of null channels has also been emphasized in 
the context of searches for un-modeled gravitational wave 
signals in ground and space based detectors. For exam- 
ple, in the LIGO- Virgo searches for un-modeled gravita- 
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tional wave bursts, it has been shown that the sensitivity 
can be improved by using the sums and differences of the 
output from the two (H 2) and four (HI) kilometer detec- 
tors at the Hanford site. The null channel H- = H1—H2 
is insensitive to gravitational waves, and has proven use- 
ful as a tool to distinguish between instrumental artifacts 
and gravitational wave signals [23|, [24[ . 

It may therefore seem a little surprising that the null 
"T" channel plays no privileged role in the present anal- 
ysis. The reason is simple: when using a Gaussian like- 
lihood function the coordinate transformation in signal 
space that produces the null channel leaves the likeli- 
hood unchanged. For example, when we repeat our anal- 
ysis using the cross spectra for the {X, Y, Z} channels we 
get results that are identical to what we found with the 
{A, E, T} channels. It is only when the instrument noise 
is not well understood, and there are significant depar- 
tures from stationarity and Gaussianity that null chan- 
nels become important. It would be naive to assume that 
the LISA data will be perfectly stationary and Gaussian, 
and we expect the T-channcl will play a key role in detec- 
tor characterization studies. That is our motivation for 
the calculation described in the next section, where we 
derive a new version of the T-channel that is insensitive 
to gravitational waves for un-equal arm-lengths. 

Bayes Factor For A,E, and T Channels vs. X,Y, and Z Channels 
100000 r . . . . . . . — . — i j 

A,E,T Channels 

30 

X Channel 

10000 - Y Channel : ! 

Z Channel : , 

1000 - / ' i 




0.01 I . . . . . . . . — I 1 

1e-13 1e-12 

"a" 

FIG. 8: Bayes factors for the X, Y, and Z channels 

With the current baseline design for LISA, the failure 
of one proof mass leaves us with a single interferometry 
channel. In most instances, it is then possible to config- 
ure the array to produce a single X type TDI channel 
(alternatively we can work with a single Beacon, Mon- 
itor or Relay channel). With one channel the number 
of cross-spectra drops from 6 to 1, and the null direc- 
tion in signal space (the T channel) is lost. On the other 
hand, the number of noise variables drops from 12 to 8 
(the effective dimension of the noise model drops from 
6 to 4 since only the sums of noise contributions along 
each arm can be inferred from the gravitational wave 
data). To study these competing effects we repeated our 
analysis using a single X type TDI channel. We per- 
formed runs on each of the X, Y and Z channels to see 



how the particular noise and signal realization in each 
channel impacted our ability to detect a stochastic back- 
ground signal. The results of this study are shown in 
Figure 8. We see that the detection threshold for a sin- 
gle channel is roughly a factor of two worse than when all 
links are operational. This result appears to contradict 
the usual statement that LISA can only detect stochas- 
tic backgrounds when the null channel is available, but it 
should be remembered that we are using the very strong 
assumption that the instrument noise is stationary and 
Gaussian, and the individual noise sources have known 
spectral shape. With these assumptions we find that the 
high frequency portion of the spectrum fixes the shot 
noise levels to very high precision, so any deviations in 
the total spectrum in the 1 — > 10 mHz range caused by a 
scale invariant stochastic background stand out in stark 
relief. 

In future studies it would be interesting to see how the 
detection limits are affected by relaxing the assumptions 
in our noise model. It would also be interesting to have 
a better sense of how well component level engineering 
models, ground testing and in-flight commissioning stud- 
ies can constrain the instrument noise model. 



V. NULL CHANNEL FOR UNEQUAL ARM 
LISA 



Unequal-Arm LISA Sensitivity Curves 





A,E Channel — 




T Channel 




New T Channel " 




il gw = 1.0*10 ra 








>v *Nj 







1e-05 0.0001 0.001 0.01 0.1 
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FIG. 9: Sensitivity curve for unequal arm LISA. Our new 
T channel restores the usual low-frequency insensitivity to a 
stochastic gravitational wave background. 

The various null channels that have been identified 
for LISA - the symmetric Sagnac channel; the Sagnac 
T channel and the Michclson T channel - arc only null if 
the arm- lengths of the detector are equal. The orbits of 
the LISA spacecraft cause the arm-lengths to vary by a 
few percent over the course of a year. In practice the arm- 
lengths will never be equal. For unequal arm-lengths we 
find that the Michclson T channel has exactly the same 
sensitivity curve as the A or E channels, at least for low 
frequencies. We found that it is possible to restore the 
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relative insensitivity of the T channel by forming a new, 
time delayed combination of the X , Y and Z channels. 
Since the arm-lengths will be approximately equal we 
write 



L a = L(l + e a ) 
L b = L(l + e b ) 
L c = L{l + e c ), 



(29) 



with \ei\ <C 1. Here an "a" subscript denotes the "12" 
arm, "b" the "23" arm, and "c" the "13" arm. In the 
frequency domain we define a modified T channel: 



T = \ (X + aY + pZ) 



(30) 



Working to leading order in and expanding a and (3 in 
a Taylor series in ///*, we are able to set the response 
function Rtt to zero out to order / 8 (the same as in the 
equal arm-length limit) with the coefficients 



jo + iji 



J2 



3-1 



k 2 4- 



(31) 



where J0J1J2 j4,&o, and fc 2 are given by: 



Jo 
ii 



(~b 



\{e a - e b ) 2 + \{e b - e c ) 2 + l(e a - e c ) 2 



32 = g(£b - e a 



17 



Ji 
fc 2 



3780 

e c - £6 
1 



(e b + e c - 2e a ) 



{(■b - e c ) ■ 



(32) 
(33) 



The a and (3 coefficients define delay operators in the 
time domain. 

Figure 9 shows the sensitivity curves for unequal arm 
LISA. We see that at low frequencies, the original T 
channel has identical sensitivity to the A and E chan- 
nels, while the new T channel restores the usual low- 
frequency insensitivity to a stochastic gravitational wave 
background. 



VI. NEXT STEPS 

We have shown that with a well constrained instru- 
ment noise model, and absent astrophysical foregrounds, 
LISA could detect a scale invariant stochastic gravita- 
tional wave background as low as fl gw = 6x 10 -13 with 
just one month of data. One line of future study is to 




5e+07 6«+n? 



FIG. 10: The time-domain ^4-channel response to the full 
galactic foreground (red); the unresolved component of the 
galactic foreground (green) and the instrument noise (blue). 
The signals were passband filtered between 0.1 — 10 mHz. 



investigate how this limit is affected by weakening the 
assumptions that went into the instrument noise model. 
Another line of future study is to sec how the limit is 
affected by astrophysical foregrounds, such as the confu- 
sion noise from un-resolved white dwarf binaries in our 
galaxy. 

It has been shown that ~ 20, 000 of the brighter galac- 
tic binaries can be individually identified and regressed 
from the LISA data [ill , but this still leaves behind a sig- 
nificant un-resolved foreground signal. It should be pos- 
sible to separate this signal from an isotropic cosmolog- 
ical stochastic background by using the time-modulated 
response to the anisotropic galactic foreground. Fig- 
ure [10] shows the time domain response to the galactic 
foreground components. Notice that at certain times of 
the year, the galactic confusion noise dips below the in- 
strument noise for almost a full month. One approach 
would be to restrict the analysis to these quiet periods, 
and achieve sensitives comparable to what we found with 
a one-month data set. A more sophisticated approach is 
to generalize our model of the cross-spectra to include 
a time-frequency representation of the response to the 
galactic foreground. The parameters for this model may 
include the bulge and disk scale heights and radii, and 
the number density of un-rcsolvcd signals as a function 
of frequency. Information from the resolved component 
could be used to place priors on this model. Work on this 
problem is currently underway, and simulated data sets 
with a full galactic foreground and an isotropic stochastic 
background have been generated for the fourth round of 
the Mock LISA Data Challenge [Hi . 

Estimates of the extra-galactic white dwarf confusion 
noise [H[ predict fl gw {f) - 3 x 10" 12 at / - 2 - 3mHz, 



9 



which should be detectable even with the galactic fore- 
ground to contend with. The extra-galactic astrophysical 
foreground will likely set the floor for detecting stochastic 
backgrounds of cosmological origin, preventing us from 
reaching the limits that the LISA observatory is theoret- 
ically capable of. 



Appendix: Cross Spectra 
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The noise cross-spectra are given by: 



AA* = - sin 2 



4 (Sfi + S12 + Si 3 + S 31 ) — 2 ( Sf 3 + Sf 2 



+5(S% 1 +S p 12 +S p 13 +S p 31 )+2(S% 3 +S p 32 ) (A.l) 



aa: = ^ 



2 I f 

sin — 



4(^2 + ^3 + ^31+^0-2(^3 



^32) 



■ cos(2u) 



2(^12 + &13 + ^23 + ^32) + 2(5*31 + S21) 



+l(s a 12 + s a 13 ) + 3(S& + ssj + + S% 2 ) 



(A.2) 



EE* S = t sin 2 (£) js* + Sf 2 + S? 3 + S* + (S* + S* ) (2 + 2 cos ) } (A.3) 



EE* = 



1G 



■ sin 



cos 2 f L 



S 23 + S 32 + S 2 i + S 31 + 2 cos 



na 1 era 1 Qa i na 
D 23 "+" J 32 "T '-'12 ' D 13 



no 1 na 
D 23 T O32 



(A.4) 



5*13 + ^32 + ^21 + 5*12 + Sfl + ^23 



(A.5) 



sin 2 ( f\ d _2cos/^n +cos 2 ( i 



16 



/* 



na _i_ oa . na 1 na 1 na r na 
°12 "r °13 "r °31 "T" °32 T" ^23 °21 



(A.6) 



AE1 



3^3 V/. 



2 cos 



lj (^13 &12 + ^31 — ^21^ 



(A.7) 



AE* 



16 



;sin 2 ( I. I 1 2cos I f ) I 5f 3 - + S 3 \ - 



3V3 \f> 



cos2 ( T ) ( ^13 — ^12 + ^23 532 ) + £31 ^21 — ^23 + ^32 



(A. 
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AT* = - si 



2f / 

/. 



1 — cos 



Sf 2 + Sf 3 + S21 + £31 — 2 ( 



'32 "T" °23, 



(A.9) 



AT* = 



■ cos(2m) 



5 a 
13 



3 



°23 



na 
°32 



3l i 



2 I §23 + "5*32 



31 "T °21 I o °12 T D 13 I o D 23 ' °3: 



_ I na 1 na \ I nil 1 na 

I D 12 T 13 I 31 -|- 2 l 



(A.10) 



et: 



(A.11) 



ST* 



16 

3^3 



sin 2 M- cos 2 



D 23 ~ D 32 ' D 21 — D 31 
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— I o 32 2 3 "r J 12 °i; 



3l i 
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